Sei ¥ + 


pcmmrnstse guetta Nentingemrincreat ESTEE 
: . » ’ ee 0 . - ~ teen Pichon ne ape ATARI NAR a eRe NS Mee oats eee ( aieaaemaniadiaa cae is 
a a a Ame a a Re eta a a ata Sanat: en TR Rt te omen we ML pews se s ; ers 
aad 1 


“ 
, N4 


. ’ wo 


y 


ae rfoe ‘ 


4 Spin eee . * ‘ 
schon » Soa, & : “ ad 
ne NT Ping PS VAG “Gaon 


So HF 


5 
Nanas Sat a nso tenn 


Si 


° re . 


Pati mn mim 


yt 
“oO, 
4 


ae 


ca 


nnn han bene nae nance maton’ 


g 


LS 


& 


Fes Ze 
Hanuts 
ke 2 


ce 


ca 


3 


This Gar 


4 


say 


b 


sey ot 
eg LS 


TR 
eas 


a nat a 
TESS Be a ee Blt. <6 eee 


Ca SS Co cm fh 


aE, {tS 


$i, 


— a Aig SABE OTs Men TC Oe 
=e eg eee ewe ee 2 7 
= . ica al 


C3 oe oo ee 


ie oer 
ise 
~~ 


" PHYSICAL AND THERMODYNAMIC 
. PROPERTIES OF POTASSIUM 


R. T. Caldwell 
D. M. Walley 
The Garrett Coevoratten 
AiResearch Manufacturing Division 


: 
Hie 
eo) 


This report was prepared by the AiResearch Mamifacturing 
Company of Arizona, a division ef The Garrett Corporation, 
Phoenix, Arizona, 28 a pert of the SHAP 50/SPUR liquid-potassium 
Rankine-cycle power-system program. The activities discussed in 
this report vere initjated uncer United States Air Force Con- 
tract AP33(615)-2289 iBPE&N: 5(6399-657)63409124). This cantract 
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This report presents an up-to-date compilation of the 
physical and thermodynamic properties of potassium. 


Thie report is assigned supplementary report number 
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ABSTRACT 


fhis report is an up-to-date compilation of potassium 
ohysical and thermodynamic properties, Data for the salid 
phase le presented, where available, over e temperature range 
frov; OOF to the mekting point. Liquid and vanor physical 
prceperties aré tabulated, in most cases, from tha melting 
point to 2400°F. The thermed compilation extends irom 
740°R to 2340°F and includes superheated vapor properties 
evaluated at pressures down te 0.2 psia. nese properties are 
presented both on a tabie an¢ cn a Hollier diagram. Listed 
properties inciude density, viseosity, surface tensicn, 
electrical resistivity, thermal conductivity, speciZic heat, 
vapor preséure, latent heats of fusion and vagorization, 
enthalpy, entropy, acoustic velecity, critical pressure and 
temperature, thermal neutron absorption and activation crose 
Sections, ionization potential, and relative magnetic 
susceptibility. 
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See Section IV D, page 16, for the symbols used in the FORTRAN 
program to compile the saturation and superheat properties 
appearing in the tables of Appendixes I and II. 


SYMBOLS AND UNITS 


Electrical resistivity, microhn-in. 


Temperature, °F 


Temperature, "R 


Thermal conductivity, Btu per hr-ft-°F 


Equilibrium specific heat, 


Density, ib per cu ft 


Absolute viscosity, lb per 


Surface tension, 1b per ft 


Enthalpy, Btu per 1b 


vidi 


Btu per 1b-°F 
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SECTION I 
INTRODUCTION 


Over the past several years, interest in the use of liquid 
metals such as potassium, sodium,rubidium, cesium, and others 
has grown to sueh an extent that the need for more accurate and 
more complete information on the physical and thermodynamic 
properties of these fluids has intensified. As a result, sev- 


eral independent laboratories, including the U. S. Naval Research 


Laboratory, Battelle Memorial Institute, and MSA Research 
Corporation, have undertaken considerable experimental work 
under Air Force and NASA sponsorship to evaluate these 
properties. 


Of all the liquid metals, potassium appears to have the 
greatest potential value as a Rankine-cycle working fluid. The 
specific volume-pressure characteristics of potassium vapor 
permit reasonable turbomachinery sizes and operating pressure 
levels over a wide range of cycle power output. Good cycle 
efficiencies can be attained with potassium at temperatures low 
enough to permit the use of relatively conventional materials 
of construction such as stainless steel. Substantially higher 
efficiencies and reduced system weight can ultimately be 
attained at the higher Loiler temperatures permitted by the use 
of refractory metals. In addition, potassium displays the 
favorable heat-transfer characteristics of all liquid metals, 
and its feasibility as a bearing lubricant has been experimen- 
tally verified. For these reasons, potassium cycles are under 
active development or consideration for use in space, marine 
applications, topping cycles in steam plants, etc.. which 
dictates the need for accurate potassium properties. 


An extensive review of the published literature relating 
to potassium has revealed inconsistencies in the property data, 
especially in the high-temperature region where experimental 


techniques are hampered by convainment problems and difficulties 


in temperature measurement. In some instances, the published 
data is considered reliable over the range of measurement, but 
is not very useful for cycle calculations encompassing lower 
and/or higher temperatures. 


The purposes of this study were to conduct an exhaustive 
survey of potassium property data, identify the most reliabie 
sources of information for each property, and compile the 
resulting data. In some cases, where equal weight was assigned 


to scattered data points, a least squares curve fit was utilized 
to estaolish the working equations used for the compilation. in 


all cuses, selections were made to insure that this document is 
internally consistent. 
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z fhe scope of this report is limited to the physical, thermo- 
2and Some of the miciear properties of potassium. For 
the convenience of the user, 211 charts, tables, and graphs are 


> presented in appendixes in the final section of the report. 


Preceding these a xes a Giscassion of the physical proper- 
ties and thermodyareic treatment cf potessium is given, along 
with 2 list of the references that were actually used to deter- 
mine Gr suppert the preperty dat: selected. The List of refer- 
eaces is not 2 cozplete ceapilation of everything that was 
examined during the course of the study. 


The dsta is presented in units that sre in common usage in 
the engineering com=amity of the United Stater, 


This document is intended to be used as a replacement for 
that portion of Reference 1 which ccvers potassium. Reference 1 
has fiiled a vital requirement over the past few years in stan- 
dardizing property data in order that everyone workizxg in the 
field would be using the same basic data. it is considered to 
have been very well done. However, since its publication, a 
great deat. of additional work has been done in establishing the 
property data of potassium (and other liquid metals}; according- 
iy, an updating is very mick in order. Weatherford, in Reference 
2. aiso makes this observation. Since References 1 and 2, like 
this document, do not report any original measurement or property 
data, they are not listed as a primary source of information but 
rather the original source was quoted. when applicable. In some 
cases, the data does indeed agree with Reference 1, but in most 
cases, additiorai later measurements have been made which super- 
sede this work. Reference 2 contains some updating and the sug- 
gested values in that reference for many properties sgree with 
the chcices made herein, whereas in other instances later work 
again superseded them. 


This document is also more complete for potassium than 
Reference 1 in that more properties over a wider temperature 
range are included. In particular, the thermcdynamic properties 
are included ip the form of superheat ‘ables, saturation tables, 
equations that can be incorporated directly into a computer pro- 
gram, 2nd a Mollier diagram of substantial size and accuracy. 
All other properties are reported in the form of graphs, tables, 
and equations which can be incorporated in @ computer program. 


Additional property measurements in all instances would, of 


course, increase the reliability and accuracy of all the property |. 


data reported, but several properties in particular nee. more 
work above a high confidence can be placed in the selected data. 
These include critical point values; the electrical resistivity, 
thermal conductivity, and specific heat of the solid; the spe- 
cific heat of the liquid at higher temperatures; and the vis- 
cosity and thermal conductivity of the saturated vapor. 
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SECTION il 
POINT VALUE PROPERTIES 


> C23 tr 


The sources of the point values selected are indicated on 
Table I, along with the values themselves. The succeeding 
paragraphs supply amplifying comments and reasons Jor the 
selections when appropriate--a practice that will be foliowed 
throughout. 


The melting point, heat of fusion, atomic weight, density 
o£ the solid, ionization potential, and relative magnetic 
susceptibility are essentially handbook vaiues which vere 
taken from the sources indicated and which are more or less 
pet aaa: and/or substantiated by later work to the accuracy 

ncicated. 


= 


The boiling point and heat of vaporization were taken from 
Reference 3 for reasons discussed more fully in Section IV. 
These point values are repeated here f@& convenience and for 
internal consistency. 


Cc 


The critical values are taken from the work of Grosse as 
reported in Reference 4 and other publications by the same 
author. Novikov. Reference 5, and others have estimated 
the critical point values of potassium with widely varying 
results, but the work of Grosse appeare to be better 
substantiated by both the experimental data and ty theory. 


The nuclear data in Table II is of a different nature 
than most of the data in this document, but is included 
because cf the widespread interest in the use of potassium 
with nuclear power sources. References 6 and 7 are the most 
recent sovrces known for this information. 
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SECTIO“, itl 
PHYSICAL PROPERTY DISCUSSZON 


A. Properties ef Sctid Potassium 


1. Electrical Resistivity 


&ithough the electrical resistivity of sciid potassiua 
should be a reletively easy property ta measure over a range of 
tempgeraturea, only a Limited number of measurements were found 
in the literature, and these were limited to a restricted ten-~ 
‘perature range or were single-point measurerernts. Sources of 
data used included References 2%, 9, 10, and 11. A eeretene 
line equaticn is considered sufficient to represent the limited 
dgte available. The equation derived is: 


p = 2.06 + 0.0096 t 
where p = electrical resistivity, microhm-in. 
t = temperature, °F 


This equation is eongidered valid over the range from §°? 
te the melting point and was used to prepare Table IV and 
bod? —_—a 
raguze 2, 


2. Thermai Conductivity 


Very few directly measured values of thermal conductivity 
for solid potassium were found in the literature. Accordingly, 
it was decided to estimate this property from the electricai 
resistivity by means of the Wiedemann-Franz-Lorenz relationshiz: 


Lo = ps 


where Lo = Lorenz constant 


p = electrical resistivity 
k = thermal conductivity 
T = absolute temperature 


The Lorenz constant ios liquid potassium is established in 
References 9 and 12 as approximately 2.15 watt-ctms per {°K}, 
which is approximately 13 percent below the theoretical value. 
It cannot be etated with certainty tnat the Lorenz constants for 
liquid and solid potassium are identical; also, the electricai 
resistivity of the solid is itself not weil established. It was 
considered sufficient to estimate a Lorenz constant for solid 
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potassius by using the few values of the thermal comductivit; 


(“> &S3 Gh 


of the solid reported in Reference 9 and the curve of the sulid- t 
phese electricai resistivity herein. The Lorenz constant thes é 
derived is, in the units explored in this report, 
0.252 SER ee ree ror about 18 percent below the theoretical 

-Ft } = 
value: 


The resuicting equation for the thermal conductivity is: 


ae 259.7 =t 
8.17 + 0.0Bit 


where k = therzal conductivity, Btu per hr-ft-°F 
t = temperature, °F 


This equation was used to generate the date for Table I¥ 
and Figure 2. 
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3. Specific Heat 


The primary refercuce selected for the specific heat of the 
solid was the work of EMI reported in Reference 13. The basic 
measurement was that of enthainy, from which an equatiocz of the 
form Hd = At + Bt@ was derived. For practical purposes, the 
specific heat is the first derivative of this equation, and the 
resulting equation in the waits used herein is 


C, = 0.1631 + 2.328 ¥ 1o *t 


ae 


E 
tr 
ip 


where C Specific heat, Btu per 1b-°F 


- 
u 


temperature, °F 
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This equa“ion was used to establish the values in Tabie IV 
ang Figure 3. The results reported in Reference 18 agree 
closely at low temperature, but have a slightly higher slope, 
with the resuit that they are approximately 1.7 percent higher 
as the melting point is appreached. The difference is not con- 
Sidered significant. 
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B. Properties of Liquid Potassium 


i. Density 


ices es on a 


or 


; Density measurements have been made by many investigators, 

4 and good agreement exist2 among the various sources of data. 
The follc.ing density equation from Reference 3 was selected 

as the recomws.nce : 
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Sthough there is some scatter, the ressite st otter inves- 
stots bracket this curve fairly weil cwer the temperature 
range from 206°F te 2000°7. Zn particeler, che resuits sf 
_ Eovikov (Reference 5} are somewhat higher, wheress these of 
- Reference 23 and the older omes of Hererence i tre icver. 


Thercal condectivities sere caicnlated aa a check from the 
@isctrical resistivity reror in the vprevicus section tbr means 
of the Siecdemom-Preoe-Lerezz reizticnship. ‘The i 
valine o£ the Lorem constant az reported by Grosse (Zeference 
271, which closely ggrees with thet of Zefezence 20, was used. 
These calculated values of thers? conductivity egree closely 
with values obtained from tne above equation. 


& Specific Fear 


Humerous imvestigetors have seée celorizetric seasurezents 
on liguid potassixa in oréer to ascertain the enthalpy and, in- 
Girectiy, the saturated licuid specific heat es fmections of 
temperature. Eelox 1000°F, there is gocd egreement between 
earlier work at the Hationel Bureas cf Standards (Reference 14} 
end more recent seasuresents at Betrelie Mexorial institute 
{Reference 13). Both sets of data can te correlated by para- 
bolic equations; hovsever, there is ¢ pronewunced deviation be- 
tgcen the twa cerves at higher texperatures, with the experizenteal 
E#i values exceeding the extrapclated HBS curve in magnitude. 


Three isolated data points obtained at the Hevai Research 
Laboratory {Reference 3; tend to confirm the exirapoiated BBE 
Gata up to 2000°F. However, the HEL investigators noted that 
there was a probable exror of 23 percent in their results, due 
to nenreproducibility in calorimetric measurements. The low 
precision is believed to stem mainly from the lowe sempie-to- 
container thermal capacity ratic. 


The need fur thick-walled cuntainers, coupled with the 
probiem of container oxidation, has rendered precise determina- 
tion of liquid-heat content extremely difficult at elevated ten- 
peratures. For this reason it is considered inadvisabie toa 
compte saturated and superheated vapor properties via the liquid 
path. For the purpose of < scing liquid-potassiu: specific 
heat up tc approximacely 1S0C°F, with cautious extrapolation to 
higner temperatures, the foliowing equation based on Reference 14 
data is recomended: 


C, = 0.2023 - 4.33 x 16°? t-$2.3 2 10°C" 
where Ce = specific heat, Btu per lb-“F 
t = temperature, °F 
Specific heat ie preserted as @ function of terperatuve 
froc. the mel-ing point to 2400°F in Table +: and rigure G3. 
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Cc. Properties of Seterated Fotassius Varor 


1. Vapor Pressere ; 
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2fter en extensive review of potassium vepor pressure in- 


¢ formstion reported in the literature. the equation Fitt the 4 
He smoothed E2L deta (Reference 3) was chesen for this compilaticn. ig 
The recommended equation is of the three-term Kirchhoff type ae 

{log Poor = 2@t b/T +c log 7). Aiter substituting the americal 5 

i constants and solving explicity for ¥5.., the equation becozes: 3 
fee rg veer £ 

nm 1.0715 x 19” i 3 
fj “Sat me 15. 58292 e ¢ 
{7} me Fi 

where Po,,- = vepor pressure, psiz Fe 

- ‘ 3 

sss f= absolute temperature, 2 


~ 


Although the equation wes primarily intended to correlate 
experinental data above the atmospheric toiling point (1394°F), 
it is believed to be sufiiciently accurate for extrapolation 
down to approzinereiy 800°F. The vapor pressure curve recon- 
mended by BMI (Reference 2G} crceses the RRL curve et epprozi- 
ms ly 10G3% ” tai ~ lu the =? 8 1 “97° 
mately 10G>°P, yielding values t are about & percent lower 
at S00°F anc 5 percent higher et 2200°F. Since the RRL therm- 


dynamic proper=y? data for notearesiim had heen selected for usa 
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elsewhere in this report, it vas deemed advisable to utilize 
the RRL vapor pressure equation for internai consistency. 
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fhe recommended vapor pressure deta is depicted as a func- 
tion cf temperature (expressed in degrees Fahrenheit} in Table 
Vili and Figure ic. 
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Specific Volume 


} The saturated vapor specific velume was computed .rom the 
4 awe Py 
equation: 


f 4 
0.2 REG T . a oe 2a s| z 
Li. Vg = PL Et 8 Pour t 3 (Paggl® + 8 (Pgag) . 
SAT fe e 
3 € 
{i where Yo = saturated vapor specific volume, cu ft per 1b i 
Se t 
a. Peat = vapor pressure, psia ; 
if Tf = absolute temperature, °R : é 

8 = second virial coefficient, psia * 
, = third viriai coefficient, psia * i 
: £ = fourth <irial coeffictent, psia > : 
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The above equation will be reccgnized as a virial equsticn 


-of stete. Yee virisi coeificients 2, y, and 6 are functions of 


pai sar et conpatibie Sith the P-¥-T deta of Reference 3. The 
ocerming tions for the viri2i coefficients are given in 
Section i¥-D of this report. 


Yapor specific voluse is listed in Table VIZI from 8C0°F to 
35 pictted Versus temperctere in Figure 11. The spe- 
cific volume of the superheated vapor can be obtained from the 
computer printout sheets in Appendix ff. 
3. weer af Varerication 
The heat of vaporization was caiculated from the Clapeyron 


equation, using the HRL vapor pressure-temperature relaticnzhip. 
Tne resulting expression is: 


Bee = + [ssarpr.ce a 0.53055 | (Ve = Yo) 


where Eee = heat @E vaporization, Btu per 1b 


T = abolute temperature, °R 

Posy = Vapor pressure, psia ~ 
Ve = saturated vapor specific volume, cu ft per lb 
Vp = saturated liquid specific volume, cu ft per lb 


The saturated liquid specific volume is simply the inverse 
of the liquid density discussed eariier. A comparison of the 
heat of vaporization calculated by different investigators 
reveals rather wice discrepancies at the high- eompere ture end. 
For example, the BMI curve (Reference 20) is 18.5 percent Lower 
than the NRL curve at 2200°F, due primarily to aif toxences in 
the experimental P-V-T data. 


Heat of vaporization computed from the above equation is 
presented in Table VIII and Pigure 12. 


4, Viscosity 


No experimental vapor viscosity data for sotassium was 
located by thie authors at the time this report was written. 
Hence, it was necessary to rely on tha foliowing theoretical 
equation suggested by Grosse (Reference 22) for the purpose of 
this compilation: 
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3 = 0.001577 /T 


Where y = viscosity, lb per ft-hr 


T = absolute temperature, °R 

The recomended equation is based on kinetic theory and is 
strictly applicable to the low-pressure monatomic gas only. 
However, to a first approximaticn the presence of small concen-~ 
trations of dimer and tetramer molecules should not appreciably 
change the viscosity. Weatherford (Reference 1) employed a 
Simpilar type of equation ts calculate the viscosity of saturated 
potassium vapor. Hewever, his values are approximately 32 per- 
cent lower, due apparently to the use of a larger-ztomc colli- 
Sion diameter. — 

Vapor viscosities calculated from the Reference 22 equation 
are listed in Table YIEIL and plotted in Figure 13. 


~ 


5. Thermal Conductivity 


It was also necessary to calculate thermal condwctivity 
data for the saturated vapor, since adequate experimental data 
was lacking. The tsermai conductivity of the monatomic gas can 
be readily estimated, like the viseczity. from kinetic theory. 
However, the saturated vapor consists of a mixrure of chemically 
reacting molecular species over the temperature range of interest, 
with up to 12 mol percent dimer and tetramer content. Hence, 
the thermal transport mechanism involves more than simple trans- 
lational energy exchange between molecules. Since an appreciable 
fraction of the heat transfer may be due to dissociation, reccm- 
bination, and diffusion of molecules through the boundary layer, 
the "equilibrium" thermal conductivity is considerably higher 
than the "frozen" thermal conductivity. 


For the purpose of this compilation, the equilibrium thermai 
conductivity was estimated from the equilibrium spe.ific heat, 
by use of the following equation: 


where k = thermal conductivity, Btu per hr-ft~-°F 


viscosity, 1b per ft-hr 


M = vapor molecular weight, 1b per lb-mol 
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equilibrium vapor specific heat, Btu per lb-°F 
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; The above equation is based on the work of Zucken and Maz- 
- well for estimation of the Prandtl number of a pSiyatomis gas. 
Fhe specific heat and molecular veight were obtained from 
Reference 3 and the computer data described later. 


Calculated thermal conductivity for the saturated vapor is 
_ Presented in Teble VIII and Figure 14. 
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SECTION IV 
THERMODYNAMIC PROPERTY DISCUSSION 
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& Review of Publisi:d Bata 


Several attempts have been made recently to establish the 
thermodynamic properties of potassium liquid and vapor over an 
extended range of temperature and pressure. Properties of 
interest include the specific volume, enthalpy, entropy, and 
specific heat of the vapor, plus a number of comparable proper- 
tiee for the liquid phase. Unfortunately most of these efforts 
have met with only limited success because of the lack of 
reliable experimental data, especially at high temperature. 
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One of the earliest published compilations of potassium 
thermodynamic properties (Reference 1) treated the vapor as an 
ideal gas mixture of monomer and dimer species whose composition 
varied in a predictable manner with temperature and pressure. 
Saturated liquid enthalpy at elevated temperature was estimated 
by extrapolating Bureau of Standards experimental data 
(Reference 14), and the saturated vapor enthalpy wae calculated 
at each temperature by adding the heat of vaporization to the 
liquid enthalpy. Internal consistency among the various proper- 
ties was secured by adjusting the constants in the vapor pressure 
equation, 


: 
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The high-temperature data of Reference 1 was in error pri- 
marily because the monomer-dimer model could not adequately 
account for the additional degree of polymerization taking place 
at high vapor pressures. In recognition of the need for more 
accurate potassium-vapor properties, several independent projects 
were initiated to measure the specific volume of potassium vapor 
at elevated temperatures along with the vapor pressure, density, 
and heat content of the saturated liquid. Unfortunately, the 

eneral lack of agreement in the experimental resulcs thus far ; 
as oe it necessary to scrutinize all pubiished data very L 
carefully. 
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After 4 cGetailed survey of experimentally determined high- 
temperature potassium thermodynamic properties, it was concluded 
that the Naval Research Laboratory data (Reference 3) is the 
best information currently available covering the temperature } 
range from 1400°F to 2500°F. Due to the absence of reliable 
high-temperature liquid enthalpy Jata, the NRL investigators k 
wisely chose to compute enthalpies and entrcpies via the monomer | 
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gas path. The P-V-T data for potassium vapor was correlated by 
a virial equation of state of the following type, which treats 
the vaper as an imperfect monatomic gas; 
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The existence of several molecular species in the vapor com- 


position ig generally accepted, The thermodynanic treatment of 


petsesium at 2 ego-called ideal gas mixture vith shifting chemical 
composition is perfectly satisfactory provided that 211 molecular 
apeciés can be identified and the entropy change associated with 
‘Bhifting composition is taken into account. The thermodynamic 


' gohfitetions made for this report were based on a virial equation 


of atate dérived froa actual P-¥-T measurements, which treats the 
veper ag an impérfect gas having constant molecular weight. All 
vapor properties were compiled from an ideal monomer gas reference 
‘state, In effect, this method avtomatically tekes into account the 
situation that the potassium is in fact neither a monomeric nor 


‘an ideal gas. he virial method was favored in the present case 


because it expedited the mathematical computations and did not 
require precise knowledge of the molecular species present in the 
vapor. Other investigaters, including References 3 and 20, have 
taken a similar approach and have demonstrated essential agree- 
ment between the virial and quasi-chemical methods at the 
preseures within the realm of interest. it is not necessary to 
compute an entropy of mixing as an explicit step when calculating 
properties by means of virial coefficients. In effect, the entropy 
change is automatically taken into account when cerrecting for gas 
imperfections. Idkewise, the vapor specific heat accounts for 
energy changes associated with shifting composition. 


B. Thermodynamic Computations Employed in This Study 


Tne RRL equation of state requires an iterative solutioa for 
specific volume. Zn order to eliminate the need for an iterative 
computation and simplify the ensuing equations for enthalpy and 
entropy corrections from the ideal gas state, the NRL data was 
refitted to a virial equation of the form: 


ep = 1+ BP + yP? + 6FS 


Theoretically, more virial coefficients should be required 
with the latter type of equation to accommodate the data. However, 
judicious selection of the preesure exponents enabled the above 
equation of state to correlate the NKL experimental data with 
sufficient accuracy. it is noted that a precedent exists in that 
& similar form of equation with one additional virial coefficient 
was employed to compile the steam tables. 


The coefficients y and 5 are tempsrature-dependent functions 
of significance only at saturated vapor temperatures above 1400°F. 
The governing equations for -y and 5 were therefere obtained by a 
curve fitting the tabulated NRL data. 


The second virial coefficient, 8, was expressed as a power 
series expression in reciprocal temperature. A leav.-squares 
curve fic was employed to evaluate the numerical constants in the 
equation. Input values of © at the lower temperatures were 
estimated Ly applying the NRL vapur pressure equation in con- 
junctien with the mcnorer-dimer equilibrium constant data from 
Reference 20. Li 
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it is noted that the NRL vapor pressure equation has been 

extrapolated several hundred degrees below the temperature range 
of the experimental data upon which it is based. for the purpose 
of computing saturated vapor compressibility and heat of vapori- 
zation in the low-temperature region. The use of a singie vapor 
pressure equation over the entire temperature range is desirable 
from the standpoint of internal consistency and continuity. 
Moreover, the NRL curve appears to adequately describe the 
scattered experimental data within the range of probable experi- 
mental error between 800°F and 1400°F. 


It is believed that the saturated and superheated vapor 
properties of interest to the turbine designer can be computed 
more accurately via the monomer gas path. However, internal 
consistency throughout a wide temperature range is of importance 
for cycle calculations involving heat balances. The predicted 
enthalpy-temperature relationship for the saturated liquid is 
considered to be a good test of internal consistency for the 
equation of state and vapor pressure equation when computing 
properties from the monomer gas base. 


There is good agreement in the liguid specific heat data 
reported in References 13 and 14 up to about 1000°F. This data, 
in turn, is closely matched by the output data from the computer 
program developed for this study. At the present time, there ic 
no way of making a precise check on internal consisteney of the 
thermodynamic equations at higher temperature by the technique 
descrited above, due to the lack of reliable calorimetric data. 


C. Description of Computer Program 


A FORTRAN program has been written to compute and sia the 
thermodynamic properties needed for cycle calculaticnas. The 
saturated liquid and vapor properties are tabulated first 
‘Appendix 1), followed by the superheated gas properties 
(Appendix ii), These tables are the actual computer srint-outs 
and are included in this report along with a Moliier chart con- 
structed from the computer data. 


The starting equations in the computer program include 
those for the vapor pressure, saturated liquid specific volume, 
and ideal monomer gas properties taken from Reference 3. These 
are followed by AiResearch-derived equations for the virial co- 
efficients and their first and second derivatives. Finally, 
the real-gas and saturated liquid properties are computed from 
the ideal monomer gas base. 


ie taneeagensa ts 
py are 


= ay Cetngyeene -- 


OO TIE TAEDA 


pan a 


SarneNpeR aewarercy 


trae aes 


eee one 


ane 
sie 


AAPL? 9 Rot EG Nn Sed eet SRDS SON Ry LOE REE EN ERR EGA FG BORE HB Uied FEASY 
2 a : a tt Sarees g ay ast = a i 
NER as EAR OR oa toas Na ieee Cle cath 


% . ANE: 


Re 1) Neh SSRIMRE mabe AAP be EIR BE Tt, PENETRATE TPN TINE a EL OOET 


Ds 


‘ 
| eer ear commer eee USeRNRNRNN § LEER TTT LEN {OTIS HR TUE) EA LRAT A NOREEN A TITEL Hi REE | WARNER DROUIN | RE MTUUARNIROON Lapel rephaNR S deme teaming parted by teneUNe eng t peRRRETELENR 


Be Computer Eousticns 


f'pzee c& Svabols | 


i : 


cd 


Superheated vapor 
Supcrheated vapor 
Superheated vapor 
Superheated vapor 


Superheated vapor 


Taput Quantities 
Temperature 


| Pressure 
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Sare,gtica pressure 


Saturéted vapor entrepy 
Enthaloy of vaporization 
Entrosy; ci vaporization 
Saturat2d liquid exthalsy 
Saturated liquid entropy 


specific volume 


enthalpy 
entropy 
specific 


acoustz~ velocity 


Safurared liquid specific volume 


Second virial ceefficieat 
third virial coefzicient 
Fourth virial ceefficient 
Sa>zrated vapor specific yo 
Satevated vapor enthalpy 
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Saturation pressure: 
7 18,727.22 
1.9715 x 10% | - fe.TET.88 
Psat = rate Sszee © 


Ideal gas enthalpy: 

~ 39.375 
H° = = 698.95 + 0.127T + 24,836 e T 
Ideal gas entropy at 1 atmosphere pressure: 


31, 126 


- ———— 


S° = 0, 18075 + 0.29283 logioT + 0.7617 € T 


Ideal gas specific heat: 


28,070 
co° ~ 0.1270 + 2.888 e ~ 


Saturated liquid specific volume: 


i 


Ve = 


Second virial coefficient: 


‘ 4 ; 5 107 
_ = - 0,8082855816 + 1336763701 %.10*)~ 9 28118 848 x 


3.885477962 x 10** 9. 747784721 x 10** | 1.548462407 x 10°° 
Bj Ts = T4 Ts 


56.099 - 6.9828 (hag) - 0.5942 (zhaq)? + 0.0408 (qhe5)° 


(1) 


(3) 


(4) 


_ 1,520716767 x 107? , 8.441706354 x 1022 _ 2.033169468 x 107° (6) 
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First cerivative of second viriel cosFficient: 
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ss 8 1.338763761 z ict | 1.9163317¢ x 20° 1. 165635289 x i5** 
a ~ @z e TT = Ts 
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, B-899112886 x 1G*5 _ 7.7%2332035 x 10°* _ 9.1ek300Ke x 102* 
Ts ft ° Pa 


_ 5.909103225 x 102 _ 7. 626551574 x 1027 
va) ‘ ee) 


ats 
Second dexivative of seccad virial coefficient: 
d®8 _ 2673527402 x 10*  5.748699509 x 10° _ &.662575554 x 10*= 
aie Ts 7 Ts : Ts 
_ L.Gkoaségke x 10°** | 4.645387221 x 107®  6.387010421 = 1672 
Ts . T7 Ts 
, E72T355558 x 1075 _2.463896417 x 102° se 
Wig:) 7 TiO Vo; 


Third virial coefficient: 


Let y ~ T-1840 


y - 2.791 x 10 “ye 
(9; 


First derivative of third virial coeffi :ient: 


d _. 7(4.5737 x 10°9y + 4.3404 x 10°7y? - 7.6743 x 10° *°y?) 
SH = 2.791 x 10°78 e 


x E ~ 4.5737 x 10 Sy - 8.6808 x lo 7y” + 2.30229 » 10°°y2 | (10) 
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Superheated vapor specific heat at constant pressure: 
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--, - SMMBURY OF POTASSHDM POIRE VALUE PROFERTIES 


185 .8°F 


25.7 Btu per ib 


Heat of fusicn 


Stouie weight 3$.100 31 
Density of the Sotid 
at SB°r- 53.7 ibs per cu ft il 


Ascospheric boiling point 1395 °F 


Heat of vaperization 
2t the boiling point 


Critical pressure 
Critical teaperature 
Critica. density 
Tonization potential 


*Relative magnetic 
susceptibility at 64.4°F 


&29.1 Btu per lb 
3380 psia 

3950°F 

10.6 lbs per en £t 
4.339 volts 


5.75 x 10-° 
{MKS dimensionless) 


*Defined as ee - 1 
Vv 


Where » = permeability of potassium, henry per meter 


b., = permeability of free space = 1.257 x 10°° henry per meter 
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SPECIFIC HEAT OF SOLID POTASSIUM 


FIGURE 3. 
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